Resonance fluorescence in a waveguide geometry 
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We show how to calculate the first- and second-order statistics of the scattered fields for an arbitrary intensity 
coherent state light field interacting with a two-level system in a waveguide geometry. Specifically, we calculate 
the resonance fluorescence from the qubit, using input-output formalism. We derive the transmission and 
reflection coefficients, and illustrate the bunching and anti-bunching of light that is scattered in the forward and 
backward directions, respectively. Our results agree with previous calculations on one- and two-photon scattering 
as well as those that are based on the master equation approach. 



I. INTRODUCTION 

Interaction between an atom and a laser beam tuned close to 
one of the atomic resonances leads to light emission from the 
atom — i.e. resonance fluorescence — with a rich set of spectral 
and temporal properties. Spectrally, as the laser intensity is 
increased the emitted light will develop symmetric side lobes 
around the central excitation frequency and the resulting spec- 
tral shape is called the Mollow triplet [1]. Temporally, the 
light emitted will also show anti-bunching with a second order 
correlation that has a minimum for zero time delay [2] . 

Recent advances in integrated optics [3, 4] and supercon- 
ducting circuits [5] make it possible to think about quantum 
systems connected to each other via waveguides that operate at 
optical or microwave frequencies. For such waveguide embed- 
ded systems, the Mollow triplet was observed in the emission 
spectra from a single superconducting qubit [6] and correlation 
measurements were also reported [7-9]. These structures were 
later shown to work as a switch [10] or a router [11]. In the op- 
tical domain, resonance fluorescence was modeled in photonic 
bandgap waveguides [12] and experimentally investigated in a 
system where a fiber was coupled to a quantum dot [13, 14]. 

Conventional modeling of resonance fluorescence focuses 
on light that is emitted in a direction perpendicular to the 
direction of the laser excitation [15-18] which results in anti- 
bunched statistics. In multi-qubit systems it is possible to ob- 
serve both bunching and anti-bunching due to the interference 
of light emission from different qubits [19-21]. In a waveguide 
geometry excitation and observation directions are co-linear as 
shown in Fig 1 . The transmitted amplitudes have contributions 
from both the incident waves and the emitted waves from the 
atom. The resonance fluorescence effect is therefore different 
from the conventional situation. In a previous study based on 
two-photon scattering off of a qubit embedded in a waveguide, 
bunching and anti-bunching of light due to the interference 
of the incoming light with the scattered fields in the transmit- 
ted and reflected directions, respectively, was predicted [22]. 
In this work, as an original contribution, we will extend the 




FIG. 1. (Color online) Schematic of the system under investigation. 
A right-going coherent state {r\^) at frequency k and an arbitrary 
intensity, propagating in a waveguide denoted by the long horizontal 
lines, is incident on a two-level system with energy separation Q. and 
a spontaneous emission rate r"'. After interacting with the qubit, 
the transmitted {r^ud and the reflected (fout) light has both a coherent 
(oj = k) and an incoherent (w k) component. 



two-photon analysis to the case where the excitation is made 
with an arbitrary intensity coherent state. We will make use of 
input-output formalism [23] recently generalized to waveguide 
structures [24] to derive analytical expressions for the second 
order correlation functions for the reflected and transmitted 
fields. We will further show that the low excitation limit of the 
coherent state solutions agrees very well with the two-photon 
results. The Mollow triplet will naturally emerge in our analy- 
sis. A distinct feature of our analysis is that we can calculate 
the multi-time correlations without specifically referring to the 
quantum regression theorem. 

The outline of this manuscript is as follows. In Section II 
we will provide the necessary definitions and derive the single 
and double time correlations for one-way waveguides by us- 
ing input-output formalism. In Section III we will extend the 
analysis to two-way waveguides and derive the spectra of the 
transmitted and reflected fields. Section IV will have the analy- 
sis on double time correlations for the scattered fields where 
we compare the coherent state and the two-photon results. We 
will conclude the manuscript in Section V. 



II. DERIVATION OF THE SINGLE AND DOUBLE TIME 
CORRELATIONS 
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A system consisting of a qubit interacting with photons in 
a waveguide is described by the Hamiltonian, H - + H\, 
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where [24] 
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Here, w is the atomic transition frequency, and are the 
creation and annihilation operators for photons at frequency co, 
respectively. cr_ and cr+ are the lowering and raising operators 
for the qubit, cr, - [<t+, cr_]. V denotes the coupling strength 
between the atomic states and the waveguide modes, and Vg 
is the group velocity of the propagating waveguide mode. In 
the derivation of the Hamiltonian we make the dipole and 
the rotating wave approximations, linearize the waveguide 
dispersion around the excitation frequency to obtain the group 
velocity and assume that the photons are at a frequency in the 
vicinity of the excitation wavelength so that the linearization 
can be justified [24]. 
We set 
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where and Uout are the input and output fields defined long 
before (fo -oo) and long after (t\ oo) the interaction 
between the qubit and the the photons takes place. The two 
fields are related by 
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The value of is in general complex valued. We define 
a/i = \aii\e'^. ojr = 2|q'j-|/ ■\/7rr is the Rabi frequency. 
The expectation value of an operator, O, is given as 

{0} = (atlOlat). 

In order to describe resonance fluorescence in a waveguide, 
three classes of correlation functions will be of importance: 
ones with one operator, ones with two operators at two different 
times and ones with three operators at two different times, i.e. 
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In order to calculate these expectation values, we use input- 
output equations (3) and multiply them from the left and the 
right with the necessary terms.' We then take the expectation 
values, make use of (4) and the commutator [ai„{t'), cr_(f)] - 
for t' > t [25] to arrive at the following set of difl'erential 
equations for all three classes of expectation values (n - 1, 2, 3) 
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Through the help of the Heisenberg equations of motion and 
the definitions of the input and output fields, we can write 
the following set of input-output equations for a single qubit 
system [24] 



dcr+(0 /2 t /I , 
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In this article we will be interested in the statistics of the 
scattered fields when a coherent state input is incident on the 
qubit. We define the incoming coherent state at frequency k as 
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These are called the optical Bloch equations with radiative 
damping. For different n, the inhomogeneous term hn and 
the initial conditions at f' = f are diff'erent: by - --, 

h2 = -f<o-+(f)) and = -^<cr+(f)cr_(f)). Previously, the 
same results were derived through the help of the quantum 
regression theorem [16, 18]. However, the derivation here fol- 
lows naturally within input-output formalism. In Appendix A 
we provide the derivation of the general solution to the Bloch 
equations and in Appendix B explicit solutions for all c„ are 
listed. 



III. DERIVATION OF THE FLUORESCENCE SPECTRUM 
OF THE TRANSMITTED AND THE REFLECTED LIGHT 

Up till now, our analysis did not distinguish between right- 
and left-going waves. Indeed, the Hamiltonian we wrote was 



' For instance, to get the second set, C2(f, f'). we need to multiply (3) evaluated 
at time (' by cr+(() from the left. 
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for a chiral (i.e. one-way) waveguide. For a regular two-way 
waveguide where fields propagate in both directions, the Hamil- 
tonian has separate input and output operators for right (r) and 
left {£) propagating waves [24]. The equations of motion be- 
come 
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We can decompose the right and left input/output states as 
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and as a result arrive at the Hamiltonian, H - Hq + H\, where 
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The fields a and a are even and odd combinations, respectively, 
of the right and left propagating fields. The interacting part of 
the Hamiltonian, Hi , depends on a only and the ii dependence 
is solely in the non-interacting part, Hq. Except for an addi- 
tional term in //q,^ the two-way Hamiltonian is very similar 
to the chiral Hamiltonian in (1). Hence, we will be able to 
make use of the results of the previous section in the analysis 
of two-way waveguides. To do so, we decompose a right going 
coherent state with frequency k into two separate — even and 
odd — channels [26] 
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where r' = t/2 absorbs the V2 factor. As one can see, the 
odd channel is interaction-free and thus is an eigenstate of 
flinCO - QoutCO whereas the even channel is subject to Hi. 
Nevertheless, it is the combination of both the even and odd 
channels that lead to the right- and left-going fields. The two 
channel expectation value of an operator O is defined as 

at at at at 
«0» = <-^;^|0|^;^). 
V2 V2 V2 V2 

In order to calculate the spectral properties of the trans- 
mitted fields, we need to calculate the Fourier transform of 
^^''out(^) ''out(f + ^0)) with respect to 6t (see Fig 1). By using (7) 
we can write 
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The application of (2) with t — > r' results in 
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One can show that the two channel expectation values of opera- 
tors are the same as their single channel expectation values [i.e. 
those in (5)] except for the substitution r — > t'. The derivation 
can be made by using (6), taking the relative expectation values 
and using (8) to simplify the results. Therefore, we can use the 
steady state values from Appendix B to arrive at 
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where 
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and 



R — cort' . 



In order to calculate the Fourier transform of this expression, 
we need to know {{cr+(t)cr_{t+6t)}} for 5f < as well. By using 
the identity «cr+(f + 6t)o--(t)}} = «cr+(f)cr_(f -i- 6t)}}* we can 
see that the expectation values for i5f < are related to those 
with 5f > by complex conjugation. The Laplace transform 
results in Appendix B thus allow us to calculate the Fourier 
transform as 
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where the function P is as defined in (B4). We will use the 
non-interacting case, that is, 



4t' 



for normalization. As a result, the coherent part of the corre- 
lation function, one which is proportional to 6(a> - k), will be 
given by 
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For reflected fields we need to do a similar analysis for 
(i^lutiO 4ut(f + St)}}. By using (7) and (2) we can see that 

{{{'Lit) ^out(f + St)}} = i«cr+(f)(r_(f + St)}}. 
The Fourier transform of this term is given by 
Gf\aj) = 3^,M£ljt){out(t + St)}}] 



1 1 



^R^ 



1 



n5{(jj — k) 



(W - /t) V2 H- 4 -H 

^ 7^ |f[-i(w-;t)]p J' 



We again normalize with respect to the non-interacting case, 
and obtain 
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for the coherently back scattered fields. The incoherent parts 
of the reflected and transmitted fields are equal to each other 
and are given by 
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In Fig 2 spectral features of the transmitted and reflected fields 
are plotted. Note that the results (9)-(10) agree with those in 
[6, 11]. 



IV. DERIVATION OF THE SECOND ORDER 
CORRELATION FUNCTION OF THE TRANSMITTED AND 
REFLECTED LIGHT 

Now that we have calculated various first order correlations 
and investigated spectral properties of scattered fields, we can 
start to look into the time dependent statistics of the transmitted 
and reflected fields. To do so, we will begin by investigating 
the second order correlation function in a chiral waveguide, 
gf\ given by 
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FIG. 2. (Color online) Coherent part of the transmitted (solid) and 
reflected (dashed) fluorescence spectrum for^ = ci^rt' = jO.1,2, 5) 
corresponding to the blue, red and green curves respectively. The 
reflected fluorescence for R = 5 (dashed green curve) is plotted after 
being multiplied by five. Inset shows the incoherent part (a) k case 
as depicted in Fig 1) of the spectrum with the Mollow triplet for zero 
detuning (D = 0). 



By using (2) and the results from Appendix B it can be shown 
that 



1 R^ 

<flout(0 flout(O) = --3-- 
T 



Similarly, after some algebra we arrive at the following formula 



for the Laplace transform of g, 
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Using the expression above, we can show that 
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where i/*^'^ [iA'^'] is the one-photon [two-photon] wavefunc- 
tion.^ As a result, we have shown that the second order statis- 
tics of a low intensity coherent state input and a two-photon 
input to a qubit are the same. 

For the two-mode case, we will need to calculate the corre- 
lation functions for the right-going (r) and the left-going {£) 



'See equations (43) and ( 1 20) in [22] for the one- and two-photon wave- 
functions, respectively. Note that the photons are at the same frequency, i.e. 
k = p. 
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FIG. 3. (Color online) Plots of g'^' for D = 2 and R = {0.2, 4) for the red and black curves, respectively. The dashed blue curve is the normalized 
two-photon wavefunction. (a) one-mode; (b) two-mode, transmitted; (c) two-mode, reflected case. As can be seen, the two-photon calculations 
are indistinguishable from resonance fluorescence ones for R = 0.2 but not for R = 4. 



fields as 
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After some algebra done by the help of an automatic non- 
commutative algebra system [27] we get 
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In the limit — > the second order correlation results for 
coherent state and two-photon inputs in a two-mode waveguide 
can be shown to equal each other, just like in the chiral case 
(see Fig 3). As was previously predicted, the interference 
of incoming and scattered fields leads to bunching and anti- 
bunching in the forward and backward directions, respectively. 
When R is increased, the response of the qubit gets saturated 
and there is less bunching in the forward direction but the 
reflected fields continue to show strong anti-bunching. In 
[28, 29] g^^^ was calculated for a low intensity coherent state 
interacting with a qubit in a waveguide where the qubit was 
coupled at a rate T' to non-waveguide modes as well. Our 
results supplement these previous investigations by analytically 
describing the scattering of an arbitrary intensity"* coherent 



state off of a qubit for the F' = case. 



V. CONCLUSION 

In this manuscript we used input-output formalism for 
waveguides to analyze how an arbitrary intensity coherent state 
scatters off of a qubit embedded in a waveguide. We provided 
analytical results for the spectra as well as the second-order 
coiTelation functions of the transmitted and reflected fields. 
This work supplements the previous work on two-photon cal- 
culations and shows that the two formulations agree for low 
intensity coherent state inputs. We predicted that the transmit- 
ted fields are bunched and the reflected fields are anti-bunched 
for coherent state inputs, similar to the case for two-photon 
scattering. Very recent experiments in circuit QED agree with 
these observations [30]. Functional devices — e.g. transistors 
[28], switches [10] or routers [11] — that make use of multilevel 
systems require both control signals that are in a coherent state 
basis and single- or multi-photon Fock states that carry the 
information. We demonstrated the versatility of input-output 
formalism with which one can do analysis either based on 
Fock states to calculate the full scattering matrix, or based 
on coherent states with an emphasis on coiTelation measure- 
ments. Additionally, it is possible to investigate non-linear 
effects such as the ac Stark [31] and Lamb [32] shifts using 
the methods developed to characterize qubit-coherent state in- 
teractions. Lastly, our approach paves the way to calculations 
involving higher-order correlation functions that become rel- 
evant when the qubit is strongly excited — in a recent cavity 
QED experiment asymmetry in time for g*^^ was demonstrated 
[33]. 
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Appendix A: General solution to the Bloch equations 



Substituting (A1)-(A2) in (A3) results in 



In this appendix we will provide the general solution to the 
differential equation 
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with the initial conditions given ait' - t by xi(f, f), x.2{t, f) and 
X3(f, t). Here b{t) is the inhomogeneous term, independent of t' . 
The solution method we use is the same as the one in [2, 15]. 
We begin by writing down the equations separately as 
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Once we take the Laplace transform of these equations with 
respect to the 6t variable, the convolution integrals simplify 
and we are left with 
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Integrating x\{t, t') from t' — t to t + 6t and making a change 
of variables results in 
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Using (A1)-(A2) we get 
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Likewise, for jiC2 and x-i we get 



X2(i) = X2(f, 0- 



.i^^e-'V*'X3(.-ifc) ^ . 



X2(f,f + <5f) = X2(f, f)e-(i/-'"'* - i^^e-^e-<'/-'">('+*' 



drne*'^"-' 
Jo 



(n-i))(/+m) 



X3(f, f + m). 



(A2) 



X3(f,f + (5f) = Xiit, f)e"^*' + -Mf)(l - e^^'^') 

+ iwRe-''' I dOTe7'"e'*<'+">e-7'*'xi(f,f + m) 



Jo 
Jo 



iwRe'*' I dme?"'e-'*<'+'">e-7%(f,f + m). (A3) 



These results are the general solution to the Bloch equations 
expressed in the Laplace domain. 



Appendix B: Correlation function calculations 

1. Single time correlations 

The calculation of (cr_(f)), {cr+(t)), and (cr,(f)) can be made 
by using the results from the previous section. The inhomoge- 
neous term is b(t) - - We assume that the atom is initially 
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in its ground state such that xi{0) - X2{0) = 0, xj,{0) - -1. 
The Laplace transforms of the expectation values are 
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The / — > cx) limit of these quantities is also of interest. We get 
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Here, D = (Q. - k)T stands for the normalized detuning fre- 
quency and R = wrt for the normalized Rabi frequency. 



£'S,[{o-+it)o-,{t + 6t)}] = 



1 + D2 + kR^ 
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Note that the Laplace transforms are taken with respect to 6t. 
3. Double time correlations of three operators 

We multiply input-output equations (3) evaluated at time f 
by cr+(f) from the left and cr_(f) from the right and take the 
expectation values to arrive at the double time correlations of 
three operators. The initial values are given by 

<0-+(f)cr_(f)cr_(f)) = <cr+(f)cr+(f)o-_(0) = 0, 
{cr^(t)cr,(t)crjt)} = - <tr+(f)(r_(0) = J-^Ml}l^ 
and the inhomogeneous term is 

T 2 

The expectation value of cr-(f) is at its steady state value given 
by (B6). If we compare the initial values and the inhomo- 
geneous term to the case of single time correlations, we see 
that they are exactly the same except for the scaling term 
(cr,(r) -I- l)/2 which is given by (B6). Thus, the results are just 
rescaled versions of the single time correlation ones and are 
given by 

Ls,W'"{cTAt)cT-(t + 6t)(T^m = (Bl) X (B6), 
£^,[e-'**'<cr+(f)cr+(f + 6t)cr^m = (B2) X (B6), 
Ls,[{cr^(t)cr,(t + 6t)cr4t)}] = (B3) X (B6), 

where the Laplace transforms are taken with respect to 6t. 



Appendix C: Short note on numerics 



2. Double time correlations of two operators 

In order to calculate (cr+(f)cr_(f')), (cr+(f)cr+(f')) and 
{o-+{t)o-z{t')) we use initial values att'-t under steady state 
conditions when f — » oo such that (cr+(f)cr_(f)) is given by 
(B6), <cr+(f)cr+(f)) = and <cr+(f)o-.(0) = -<o-+(f)) is given 
by (B5). The inhomogeneous term is b(t) - - ^{cr+(t)). After 
some algebra we get 



Ls,W'"{cr^(t)crJt + 6t))] = 



1+D2 + \R^ sP{s) 



_1^2g-2i(J+2iA:r 



^[..Ui(Q-fc)]-l 
s+^--i{D.-k) 



The differential equations that we analyzed so far can be 
transformed into time-independent forms by the substitution 



and &, - (Ty, 



where k is the frequency of the incoming photons. For instance, 
the single time expectation values of o"_, &+ and &^ can be 
written as 



d 
df 











( ^ ) 


<<^+> 


= M 




+ 













2 



where the matrix M is given by 



M = 



-(i+iQ) 


iwRC" 







i^e"* 



-(L-iCl) -ife- 
-iwRc'^ -I 



and = Q.- k. Other expectation values have the same form 
as well. This is a much more convenient formulation for purely 
numerical studies with which we verified the analytical results 
reported in the previous appendices. 
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